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Powerful radio-loud active galactic nuclei (AGN) with large Mpc-scale jets have been theoretically motivated as emitters of high- 
energy cosmic rays. Recent radio observations have established a populous class of young radio-loud galaxies with compact (< 1 
kpc) symmetric jets that are morphologically similar to large-scale AGNs. We show that these compact AGNs, so-called compact 
symmetric objects (CSOs), can accelerate protons up to 10 20 eV at their hot spots via a Fermi type mechanism on the assumption of 
efficient acceleration. The required magnetic field strengths are comparable to those derived from the minimum energy condition. 
We further show that the accelerated protons can escape through the photon fields of the cocoon without significant energy loss. 
However, the local number density of powerful CSOs is insufficient for CSOs to power the entire observed flux of ultra-high-energy 
cosmic rays, providing maximally only a few percent. A heavy composition of UHECRs allows more CSOs to accelerate particles 
to UHECR energies, but escaping the cocoon is difficult. We comment on a method that may test CSOs as UHECR sources. 
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1. Introduction 

Radio-loud active galactic nuclei (AGN) are extremely lu- 
minous objects displaying non-thermal activity across many 
wavebands from radio to TeV gamma rays. The radio emis- 
sion is attributed to synchrotron radiation of accelerated elec- 
trons, while the gamma-ray emission may be attributed to (i) 
inverse-Compton of accelerated electrons on either synchrotron 
photons (synchrotron self-compton, e.g., HH2]) or external pho- 
tons (external compton, e.g., [31141]), (ii) decays of neutral pions 
produced in hadronic interactions of accelerated protons or nu- 
clei, e.g., |5|,|6j] and/or (iii) synchrotron radiation of ultra-high- 
energy (UHE) protons, e.g., fl0.Hl]- Although leptonic mod- 
els provide satisfactory fits for many AGNs, the hadronic sce- 
nario is favored for some specific ones |9]. AGNs have there- 
fore been extensively studied as potential sources of ultra-high- 
energy cosmic rays (UHECRs) lid 

Most AGN gamma-ray detections are of the blazar-class of 

AGNs, which are known to have relativistic jets pointed towards 

I — I 

the observer. According to the unification model of AGNs [16], 
radio-loud Fanaroff-Riley (FR) galaxies and blazars originate 
from the same population, with the difference due to viewing 
angle. Unlike blazars, the jets in FR galaxies are not pointed 
towards the observer, and therefore the gamma-ray fluxes are 
expected to be much smaller. However, gamma rays have re- 
cently been detected from the giant lobe of a nearby FR galaxy 
ifnll . Q Considering the possibility of hadronic components in 



'Note that the first detected non-blazar source is M87 19], although the 
emission is from the core. 



blazars, hadronic acceleration could also occur at the hotspots 
and/or lobes of FR galaxies. Indeed, the hot spots of strong FR 
II galaxies have been theoretically motivated as sources of for 
some time 1 10]. 

Interestingly, recent progress on VLBI observations have re- 
vealed the existence of compact radio-loud AGNs that are mor- 
phologically similar to FR II galaxies — e.g., they contain sym- 
metric radio lobes and terminal hot spots — but the linear size, 
LS , is considerably smaller, < 1 kpc (compact symmetric ob- 
jects, CSOs). Very high radio powers are observed, up to 
L syn ~ 10 43 -10 46 erg s , and their radio power per unit fre- 
quency has a turnover at ~ GHz (GHz-peaked spectrum ob- 
jects, GPS; see, e.g., Ref. Ir20ll ). Due to these observed prop- 
erties, CSOs are thought to be precursors of radio-loud AGNs, 
i.e., observed at an early stage of jet expansion Oil 12211 . Flux- 
limited samples of radio sources find that a substantial fraction 
of the order of 10% are CSOs and related compact AGNs. Since 
this is larger than what a simple one-to-one evolution scenario 
to large-scale FR II galaxies would predict, various luminosity 
function and evolution scenarios have been proposed Ir23l 24]. 
The strong non-thermal activity, similarity to radio galaxies 
with large jets, and the significant fraction in flux-limited ra- 
dio catalogs, motivate investigating whether CSOs can generate 
UHECRs. 

In this study, we discuss hadronic particle acceleration at the 
hot spots of CSOs. The possibility of acceleration in similar 
conditions has been studied by Ref. Ir25ll . where the author dis- 
cussed the capability of accelerating protons up to 10 20 eV at 
~ 100 pc from the central black hole. We extend this by includ- 
ing the treatment of the energy source of the magnetic field, 
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co-evolution of radio jets with cocoons, and the escape of the 
protons out of the cocoon. Furthermore we put this in the con- 
text of CSOs, and discuss whether CSOs are able to explain a 
significant fraction of the observed UHECRs based on recent 
data of CSOs and UHECR observatories. We consider mainly 
protons as the composition of UHECRs, but we also discuss 
applications to heavy compositions. 

The paper is organized as follows. In Section 2 we review 
our adopted model for the hydrodynamical evolution of CSOs. 
In Section 3 we discuss our investigation of the generation of 
cosmic rays at CSO hotspots, including magnetic field energet- 
ics, cosmic ray acceleration, cosmic ray energy loss, and cosmic 
ray escape through the cocoon. In Section 4 we discuss whether 
CSOs can power the observed UHECR and discuss the impact 
of heavy composition. We conclude with a summary in Section 
5. 



2. CSO model 



AGN jets are thought to be relativistic, 



c, and are decel- 



erated by shocks produced by interactions of the jet with the 
interstellar medium of the host galaxy. In general, terminal 
shocks as well as reverse shocks (RSs) are generated, where 
large fractions of the jet kinetic energy are dissipated. The 
shocked hot plasma subsequently expands sideways to produce 
the cocoon. Substantial energy is also converted to accelerated 
particles, as shown by the observed non-thermal radiation from 
the shocked jet region near the RSs, so-called hot spots (see 
Ref. [26] for a schematic view). We consider UHECR genera- 
tion at the RSs which generates the hot spot. 

In order to physically describe the evolution of CSOs, we 
adopt the simple model worked out in Ref. [27], originally for 
classical radio-loud galaxies expanding in an ambient medium. 
This model consists of 3 equations. The first is the balance 
between the momentum flux of the jet and the ram pressure of 
the ambient medium, 



— = n a tn„v h 2 A h , 



(1) 



where Lj, n a , m p , v/,, A/, are the jet kinetic power, the number 
density of ambient medium, the proton mass, the advance ve- 
locity of the hot spot, and the cross-sectional area of the bow 
shock at the end of the cocoon, respectively. The second is 
the balance between the pressure in the cocoon, p, and the ram 
pressure, 

p = n a m p v c 2 , (2) 

where v c is the sideways expansion velocity of the cocoon (the 
speed of the shock driven by the overpressured cocoon). The 
third is energy conservation, 



Ltott 

(r-D' 



(3) 



where V, y, L tot , and t are the cocoon's volume, the adiabatic 
index of the relativistic cocoon's fluid (y = 4/3), the total lumi- 
nosity of the jet, and time, respectively. 



The Lj is related to L tot by Lj - ri^L tot , where tj^ is the 
fraction of the jet kinetic power to the total jet power. Eq. Q] 
describes the dynamics of the jet head and we assume L tot to 
be constant during the active phase of CSOs. Since we fo- 
cus on CSOs with LS < 1 kpc inside elliptical galaxies, the 
density of the hot ambient medium can be assumed to be con- 
stant, n„ ~ 0.1 cm -3 [28]. Following these assumptions and that 
l c 2 oc t Q29D, which was proposed to reproduce the initial phase 
of jet propagation according to simulations of Ref. 13011 . leads 
to Vh °c f° and A/, oc t°. l c is the transverse size of the cocoon. 
The constant v/, is in good agreement with observations, which 
typically find v/, ~ 0.1c (see, e.g., Ref. O1U32I0 . Based on this 
model, the size of a particle acceleration region is estimated as 



1/2 j 1/2 
l"-h nrv^K.0. 6^tot,46 



1/2 a 



(4) 



where ?7k,o.6 = 77k/0.6, i tot ,46 = ^tot/10 46 erg s *, n a -i = 
n fl /10 _I cirT 3 , and /?/,,- 1 — (v/,/c)/10 _1 . Note that we assume 
Ah to be comparable with the cross-sectional area of the hot 
spot. 



3. Cosmic Ray Generation 

3.1. Cosmic Ray Acceleration 

In the absence of energy-loss processes (we will discuss these 
in section [331 ), the maximum energy of particles accelerated at 
the RS can be estimated by comparing the time scale to acceler- 
ate particles and the minimum time scale in which the particles 
stay in the acceleration region [33]. The former time scale is 



9 F r' g 



lxlo io^o 



B'_ 3 r 



(5) 



where r' g , e p>20 = e p /10 2U eV, B'_ 3 = fi'/lO^ G, and V are the 
Larmor radius of the particle in the RS rest frame (RSF), the 
proton energy in the observer frame (OF), the strength of mag- 
netic field in the RSF, and the Lorentz factor of the RS in the 
OF, respectively. Throughout the paper, primes (/) explicitly 
denote quantities in the RSF. Since the velocity of the shocked 
jet in the OF is /?/, and in the RSF is 1/3 for the strong shock 
limit HS, the velocity of the RS in the OF, B RS , is 



RS 



1/3 -B h 
1- (1/3)0* 



0.25, 



(6) 



and thus T = (1 -yS RS 2 ) _I/2 ~ 1. We take F to be a con- 
stant, which is equivalent to assuming the diffusion coefficient 
for accelerated particles to be proportional to the Bohm diffu- 
sion coefficient; Of > 10 is a fairly conservative value, while 
Of ~ 1 may be achieved for mildly relativistic shocks (see, e.g., 
Ref. [36]). Note that Eq. <j3j implies that the region to accelerate 
protons up to 10 20 eV satisfies LS > 100 pc from LS /v/, > r acc 
where T acc = IY^. 

There are two time scales that can limit the time for particles 
to stay in the acceleration region, the escape time r' esc , and the 
dynamical time . The lower limit for Tg SC is given by the 
shock-crossing time in the upstream fluid, T' esc ~ r/JFB'c, while 



2 



the dynamical time scale is r' dyn ~ LS /FBrsc, where B' = 1 is 
the velocity of unshocked jet assuming the strong shock limit. 
For our focus (LS > 100 pc), Tg SC < T' dyn holds quite generally. 
Requiring T f acc < T f esc , we obtain 



(7) 



Thus, a few mG magnetic field is required to accelerated pro- 
tons up to 10 20 eV even if protons are accelerated efficiently 
(0 F ~ 1). 

In fact, magnetic fields of a few mG have been estimated for 
CSO hotspots, from the minimum energy condition (e.g., [y]\\ ). 



_ 1 67rAL hs 



2/7 



T 2/7 

= 2 — mG, 



-6/7 



(8) 



where Lhs.44 = Lh s /10 44 erg s _1 , Vh s = 47rr/, 3 /3, and A = 3.1 x 
10 7 in cgs units are the synchrotron luminosity in the hot spot, 
the volume of the hot spot, and a numerical factor, respectively. 
The factor A is dependent on the spectral index of synchrotron 
radio emission per unit frequency above the turnover frequency, 
denoted by a B37I1 . For our calculations we adopt the averaged 
value observed in CSOs, a — 0.73 [38]. The value of Eq. dHJ is 
comparable with the required value in Eq. (Q) for Op ~ \, Thus, 
the energy of protons can reach 10 2() eV if they are accelerated 
efficiently. 

3.2. Hot Spot Energetics 

The magnetic field of Eq. (O corresponds to a luminosity of 



B' 2 

L B = —4nrlT 2 B'c ~ 2 x lO 45 F V 2O ergs- 

07T ' ' 



(9) 



where dependencies on parameters of order unity — B' and T — 
have been removed for clarity. One sees the dependency on L tot 
has cancelled; in fact, Eq. © is a simple generic requirement 
of the magnetic luminosity, for a system accelerating protons to 
10 2() eV under stochastic acceleration mechanisms, i.e., those 
that must confine the accelerated particle fill l39l [Till . Since 
the magnetic field is powered by the jet energetics, one imme- 
diately sees that L tot > 10 45 5 erg s for 6f ~ 1 is required to 
accelerate protons to 10 2() eV, where Lb = T]^L iot and t]s < 1 
is a parameter. This is a large luminosity but it is not unrea- 
sonable, as we discuss next. Note also that this condition is 
relaxed for heavy nuclei, which have smaller Larmor radii for 
a given particle energy; for example, the required Lb for pure 
iron composition would be 10 42 4 erg s _1 01411 . 

The value of L tot can be derived from the radio synchrotron 
luminosity, L syn . The total energy of synchrotron radiation is 

4-o- T U B U e 4- , 

L syn = . T - r«7f/(r), (10) 
3m e c 3 

where m e and cr T are the usual electron mass and Thomson scat- 
tering cross section, Ub and U e are the hot spot magnetic field 
and accelerated electron energy densities (both in the OF), and 
y is the Lorentz factor of accelerated electrons, respectively. 
Here, 



fir) 



r^dyAy)} ' r"dyy 2 ^(y), (11) 



and dN e ldy is the spectrum of accelerated electrons. 

While shock acceleration predicts a single power-law spec- 
trum, we assume a broken power-law spectrum, 



dNe 
dy 



1 < y < y br 



(12) 



reflecting the expected suppression at hig h energies due to 
synchrotron energy losses (see, e.g., i40IP . The value of 
7br is derived by equating the energy loss time, T eS ync = 
3m e 2 cP /4o~Te e U B , to the time scale in which electrons escape 
from the hot spot, t s . The maximum y of electrons is esti- 



mated by T a , 



to be y max = 4 x lO 9 p - l/2 B^- l/2 . 



The synchrotron photon spectral slope a is related to s\ by 
a = (si - l)/2, and would steepen above the break frequency 
Vbr corresponding to y\, T . 

Although the escape mechanism is unclear at present, the 
value of t s can be limited. The minimum must be r/,/c = 



2 x 10V,o.6 I/2 £, 



*J/ 2 n , 
tot,46 "a,-l 



-1/2 



h,-\ 



yr, and the maximum is 



the age of the source, ~ LS/vf, — 3 x 10 3 L5 , 2jS/ ! ,-i _1 yr, where 
LS2 = LS /1Q 2 pc. For f s = 10 2 yr, corresponding to the min- 
imum, Tbr = 2.4 x lQrt S £~ 1 B-3~ 2 . Larger ages would result 
in smaller y\, r . Observations generally agree with this simple 
formulation. For the handful of CSOs where the hot spots and 
lobes are separately resolved spectroscopically, the Vbr in the 
lobes are lower (older) than those in the hot spots 14111 . This 
supports the framework that electrons are accelerated at the hot 
spots and subsequently expand sideways to form the lobe and 
cocoon. Note that a is observed to lie in 0.5 < a < 1 .0 with an 
average of 0.73 l38ll . consistent with theoretical expectations of 
si. 

Adopting the observed average a = 0.73 (si = 2.46) and 
Tmax = 4 x 10 9 , we can obtain f(y) ~ 10 3 in the range of 
uncertainty of t s . Now, U B = r 2 U' B and U e /T 2 U' e = 
r] e L lot /(47rrh 2 r 2 B e c), where B e and 7] e are electron escape ve- 
locity normalized by c and the fraction of relativistic electron 
luminosity to the total jet power, respectively. Although B e is 
uncertain as described above, we assume it to be the velocity of 
the downstream fluid at the RSF for the estimation. Substituting 
these into Eq. (TTOb . we obtain, 

'7B,0.2'7e,0.2itot,46 1/2 "a' ,2 1 f(y) 



^=5xl0- 3 - 

J-AsA 



/^O.A-lA'.l/3 



10 3 



(13) 



where ^0.2 = ^c/0.2, 77^,0.2 = ?7b/0.2. Thus, despite uncertain- 
ties, L syn is ~ 1% of L tot in this model. Some previous works 
find L syn /L t ot ~ 10% Il40ll42ll . which would imply a jet power 
of L tot ~ 10 45-46 erg s -1 for an observed synchrotron luminos- 
ity L syn ~ 10 44 erg s _1 . Our adopted value of the jet power, 
Ltot = 10 46 erg s~', is on the high-L end of the luminosity func- 
tion, as one would expect for generators of UHECRs. 

3.3. Cosmic Ray Energy Loss 

Here we discuss energy-loss processes that can inhibit the ac- 
celeration of protons up to the highest energies. These include 
proton synchrotron radiation and scattering with ambient pho- 
ton fields. For the photon field we include synchrotron radiation 
of electrons that are accelerated in the same region as protons, 



as well as the disk, torus, and stellar emissions. As we show 
below, protons are not prevented from being accelerated up to 
10 20 eV. 

The time scale for proton synchrotron energy-loss is 

i 4 3 

3m„ c 



4cr T m e 2 e p U' B 



2x 10 13 e 



p,20 



(14) 



adopting the magnetic field of Eq. Q. Since T acc <k Tp Syn we 
confirm proton synchrotron energy loss is not important. 

We normalize the synchrotron photon field from the observed 
non-thermal radio emission. We expect a spectral steepening at 
v br = 3 x 10 5 B_3 3 f. Sj 2~ 2 GHz JH, H corresponding to y bl , 
above which the spectral index steepens from a to a + 1/2. 
Given the observed values of a, the radiated energy peaks at 
Vbr- Furthermore, synchrotron photons are self-absorbed at 
low frequencies, and the spectral index changes to 5/2 below 
the turnover frequency v sa ~ 0.5 GHz, estimated by equating 
the optical depth of synchrotron self-absorption, ff sa r/,, with 
unity, where a sa is the absorption coefficient of synchrotron 
self-absorption given in standard textbooks, e.g., 14511 . There- 
fore, assuming all the photons are emitted from the hot spot, the 
synchrotron power per unit frequency L v is 



U = 




(v < v sa ) 

(v sa < V < v br ) 

(Vbr < V) 



where 
N 



-'syn 



(15) 



(16) 



(f-ry^ + fT^ + s-bH 

is the normalization. The average energy density of photons is 



U s 



3L, 



syn 



yn 



2x 10 



-8 



£syn,44?7K,().6£tot. 



46 



■ erg cm 3 .(17) 



We treat the photon fields due to an accretion disk, dusty 
torus around the AGN, and stellar emission in a similar way to 
Ref. [40] . Unlike many GPS quasars where strong ultraviolet 
(UV) emission have been directly observed, the UV emission is 
expected to be largely absorbed by a dusty torus in the case of 
GPS galaxies. We thus assume that GPS galaxies are similar in 
its intrinsic UV luminosity to GPS quasars, with Luv ~ 10 45 - 
10 47 erg s" 1 J4^]. The energy density of UV photons in the 
acceleration region is then 



t/uv 



3 x 10 



_7 ^UV.46 



erg cm 



-3 



(18) 



4nLS 2 c LS 2 
A large fraction of the UV energy density, parametrized by jjir, 
can be absorbed and re-radiated in the infrared (IR) by a dusty 
torus. The energy density is 



U 



IR 



„ 1n _ 8 '7lR,-0.5AjV,46 -3 

9 x 10 erg cm . 

LS 2 



(19) 



Finally, the energy density of star light photons is £7 stal = 
3Ly /Anr s 2 c, where r s is the core radius of the stellar distribu- 
tion which correlates with the V-band luminosity of the host 
galaxy L v as L v /10 45 erg s -1 ~ r s /l kpc [47]. Therefore, 



Ustar = 8 x 10 10 Lv,45 erg cm 3 . 



(20) 
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Figure 1 : Bethe-Heitler pair creation energy-loss time scale and mean free time 
of photopion production in the CSO hotspot photon field. We consider syn- 
chrotron photons, central accretion disk UV, dusty torus IR, and stellar light. 
The synchrotron cooling and acceleration times are also shown for comparison. 
The acceleration of protons is limited at ~ 10 20 eV by the time scale for parti- 
cles to stay in the acceleration (see Eq. 0. Here, we fix B = 2 mG, r/, = 70 

10 eV, L uv = 10 46 erg s _1 , <er> IR = 0.1 eV, 



pc, 



10* 



erg s 



•\ <e)uv 

n 45.5 , 



^<IR = ^IR-O.s^uv = 10 erg s , and LS = 100 pc. 



From the above estimates we see that t/uv > Ub > £Ar > 

The time scale for inverse-Compton energy loss is thus dom- 
inated by tJ cs uv ss 6 x 10 12 e / , ; 2o _1 iuv,46 _1 i52 2 s, which is still 
much larger than t^ cc . Thus, ICS on the various photon fields 
can be neglected. 

In Fig.Q]we show the energy-loss time scale of Bethe-Heitler 
pair creation and the mean free time of photopion produc- 
tion in the synchrotron photon field, the UV, and IR radiation 
fields. The spectra of the UV and IR radiation are assumed 
to be thermal distributions with the average photon energies of 
(e)uv = 10 and (e) IR = 0.1 eV, respectively. In terms of target 
photon number density, synchrotron photons play the dominant 
role, which in turn depend on the value of f s . To be cautious, 
we use the minimum f s = 2 x 10 2 yr, which corresponds to the 
direct escape of accelerated electrons. We find that photopion 
production is the dominant energy loss mechanism, but it is not 
important up to and including 10 20 eV. This is true between the 
minimum t s and the maximum t s . 

We conclude the acceleration capability is maintained under 
our adopted CSO evolution scenario; that is, until the hot spot 
dissipates or the dynamics change. The former can arise if the 
AGN activity terminates [48]. The dynamics can change if LS 
exceeds ~ 1 kpc or LS /v/, > 10 4 yr. This is because A/, is 
unchanged in our CSO model and the photon fields become 
weaker with increasing radius. When the distance of the hot 
spot exceeds ~ 1 kpc, the approximation of a constant n a is 
no long valid, and A/, is no longer constant. This phase is not 
treated in this paper because we focus on the early epochs of 
radio-loud AGNs. 
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3.4. Cosmic Ray Escape through the Cocoon 

In order to be a viable UHECR source, the protons accel- 
erated in the hot spot must survive their escape through the co- 
coon. The lobe size is characterized by the radius of the cocoon, 



lc = 



1 8L tot LS 2 
K 3mn p n a v{ : 



1/4 7-1/4 rc l/2 

~ 2xl0 2 lot ' 46 2 

- L X 1U 3 

a,-l 



pc, 



(21) 



and we compare this with various relevant distance scales of the 
proton. 

Firstly, if the magnetic energy is conserved during the CSO 
evolution and cocoon formation, the energy density of magnetic 
field in the cocoon is estimated as 



2nL 2 LS 



1x10 



-, 1/4 L, 



tot.46 



3/4 



erg cm 3 , (22) 



adopting a volume of 2nl c 2 LS . Note that the shape of the co- 
coon is approximated to a cylinder in the CSO evolution model 
described in Section [2] This yields a magnetic field strength of 
~ 3 mG field for rjs - 0.2, which is unnatural given the hot 
spot field is ~ 2 mG. In reality, the magnetic energy is likely 
dissipated or leaking from the cocoon at some non-trivial time 
scale. In addition, it is difficult to realize a uniform magnetic 
field over this size, and a strong random or turbulent field com- 
ponent is expected. If the coherent length is smaller than the 
Larmor radius of the UHECR proton, the protons propagate in 
the lobe diffusively and can escape. Even in the extreme case 
of a uniform magnetic field, one can estimate that the Larmor 
radius of a UHECR proton is larger than l c for B < 1 mG. 

The UV energy density is the largest of all energy densities 
(including photons and magnetic field; Uuv, Ujr, U^ yn , U C B ), 
and we find that ct^ s ^ 6 x 10 e p> 20 Ltjv,46 £>S2 2 pc. This is 
much larger than l c , showing radiative losses are unimportant. 
Bethe-Heitler and photopion production energy losses can be 
shown to be similarly unimportant. To be conservative, we use 
an upper limit on the number density of the synchrotron photons 
as described below. Observationally, the radio emission energy 
from the lobe (cocoon) is weaker than that from the hot spot, 
i.e., L c syn < L syn , and Vb r is smaller in the cocoon. Thus, spectral 
modeling of Eq. [B]with t s ~ LS/v h = 3 x lO^^-i" 1 yr 
gives an upper limit of the radio photons in the cocoon. Note 
that the change of v sa is small and does not affect the normal- 
ization of the synchrotron photons. For Bethe-Heitler pair cre- 
ation, the energy-loss length is at least 10 4 pc at 10 20 eV. Sim- 
ilarly, the mean free path of protons for photopion production 
in this radio field is at least 10 3 pc at 10 20 eV. Therefore, nei- 
ther Bethe-Heitler nor photopion production affect the escape 
of 10 20 eV protons. 



4. Discussion 

In this paper we focused on the early epochs of radio-loud 
AGNs, namely CSOs, and demonstrated that the physical con- 
ditions in their hot spots allow for the acceleration and escape 
of protons at ~ 10 20 eV. We now discuss the total energetics and 
the capability of confirming this scenario. 
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Figure 2: Same as Fig. [T] but for irons. The photopion mean free time has been 
replaced by the more important photo-disintegration mean free time. 



The energy budget (per volume) required to reproduce the 
observed flux of the UHECR at Earth is nominally £(e ; , > 
10 19 eV) ~ 2 x 10 44 erg Mpc~ 3 yr 1 |H EH . Uncertainties 
of the source spectrum introduce a factor 2 uncertainty to this 
value. Since UHECRs with energies above 6 x 10 19 eV travel 
maximally 200 Mpc because of interactions with cosmic mi- 
crowave background (CMB) photons ||52, 53], whether CSOs 
are the main sources of the observed UHECR depends on the 
CSO density in the local Universe. 

Since CSOs are a subclass of radio galaxies, the number 
density of CSOs cannot be larger than that of radio galaxies. 
The number density of radio galaxies with luminosities above 
L syn ~ 10 44 erg Mpc~ 3 in the local Universe is a; 9 x 10~ 9 
Mpc~ 3 15411 . The luminosity function of GPSs was estimated 
by Ref. 15511 . assuming that their cosmological evolution is sim- 
ilar to large-scale jet objects. Following this result, the number 
density of GPSs in the local Universe is estimated as 10" 10 - 
10" 11 Mpc -3 . Note that this local GPS to radio galaxy ratio is 
at least one order of magnitude smaller than observed in cat- 
alogs of high-redshift radio galaxies (23, 24]. Now, the total 
energetics of UHECR from CSOs is S — r7^Lt ot n cso , where 
fj p is the fraction of the total luminosity L tot that is given to 
protons with e p > 10 19 eV. Requiring CSOs to power the ob- 
served UHECR energetics, and that f] p < j]b - 0.2, gives that 
"cso ^ 3 x 10~ 9 Mpc -3 . The required number density is larger 
than the estimated number density of GPSs. While the uncer- 
tainties are large, adopting the local CSO number density to be 
10~ 10 Mpc 3 , maximally only a few percent of the total UHECR 
flux can be due to CSOs. 

The situation changes for heavy nuclei composition of UHE- 
CRs. A recent result of the Pierre Auger Observatory (PAO) 
indicates that UHECRs are dominated by heavy nuclei 15611 (al- 
though the High Resolution Fly's Eye (HiRes) reports a proton- 



dominated composition [57, 58]). Heavy nuclei dominated 
composition relaxes the total energy budget of UHECRs. An 
effect is to decrease Eq. dg}, and thus decrease the required CSO 
total luminosity L tot II 1 411 . The number density of CSOs above 
the threshold luminosity then increases, so that CSOs contribute 
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Figure 3: Mean free paths of photopion production of protons and photo- 
disintegration of irons in the IR photon field which is approximated as a delta 
function with erR = 0. 1 eV. 



CRs by Galactic mag netic field and intergalactic magnetic field 
(IGMF) S EE il El 0, S S (but see Refs. HI [68 ] 
for the possibility of spurious correlation). Thus, the correla- 
tion between the arrival directions of UHECRs and the posi- 
tion of CSO candidates could be a powerful way to confirm 
our scenario. Since CSOs are subdominant sources of UHE- 
CRs, much more detected UHECRs are required. Furthermore, 
the intermittency of CSO activity with the duration of 10 4 -10 5 
yr has been suggested to explain the observed overabundance 
of CSOs compared to predictions of a uniform evolution sce- 
nario to large-scale radio sources 114 811 . If the active phase is 
shorter than the time-delay of UHECRs during propagation, 
CSOs which emitted UHECRs are possibly inactive when the 
UHECRs arrive at Earth. Note that there is large uncertainty 
on IGMF and the time-delay of UHECRs strongly depends on 
IGMF modellings [Zll IZ2l llZD . In either case, we may 

not observe active CSOs, but instead CSO remnants Il73ll . to- 
wards the directions of UHECRs. 



significantly more towards the UHECR energy budget. The rel- 
ative importance of CSOs depends on the luminosity function 
of other potential UHECR sources. Although quantitative dis- 
cussions of the effects of heavy nuclei depend on the mixture of 
nuclei and are beyond the present paper, we give a brief discus- 
sion of the acceleration and escape of iron at the hotspot below. 

Figure[2]is similar to FigureQ]but for irons, showing the time 
scales for acceleration, synchrotron energy-loss, Bethe-Heitler 
pair creation energy-loss, and the mean free time for the photo- 
disintegration of iron. We find that the acceleration of iron nu- 
clei to 10 20 eV is possible. The acceleration time scale scales 
as Z _1 while the synchrotron time scale scales as A 4 Z~ 2 (see 
Equations and [T4l which work in favor for iron nuclei. In 
order to calculate the mean free time for photo-disintegration, 
we adopt the parametrized cross-section in Ref. [59]. Although 
the cross-section for iron at the first peak of the GDR is ~ 100 
times larger than the cross-section for protons at the delta reso- 
nance, the acceleration time scale is sufficiently small for iron, 
and iron can be accelerated to 10 20 eV. 

However, the 10 20 eV irons cannot escape from the cocoon. 
We demonstrate this using IR photons. Approximating the IR 
photon spectrum as a delta function at e = 0.1 eV, we calcu- 
late the mean free path of irons for photo-disintegration. This 
is shown in Figure [3] The mean free path of protons for photo- 
pion production is also shown for comparison. The mean free 
path of irons in the IR photon field is minimized at ~ 10 19 
eV due to GDR to * 10 pc. Also, this field predicts that the 
mean free path of 10 20 eV iron nuclei is less than 10 2 pc due to 
baryonic resonance. This simple investigation shows that iron 
nuclei would photo-disintegrate before escaping from the co- 
coon, of length l c =i 2 x 10 2 pc. Since photo-disintegration is 
the process of nuclear breaking, it is possible that nuclei with 
lower nuclear number and lower energy could diffusively es- 
cape from the cocoon. We conclude that while CSOs could be 
high-energy nuclei emitters, the details depend largely on the 
structure of magnetic field in the cocoon. 

The recent PAO data indicates small deflections of UHE- 



5. Summary 

In summary, we investigated CSOs as a possible source of 
UHECRs, demonstrating that their hot spots may accelerate 
protons to 10 20 eV. Although the local number density of CSOs 
capable of accelerating UHECRs contains larger uncertainties, 
it is generally insufficient to power the entire UHECR observed 
on Earth. An iron composition could increase the contribution 
from CSOs to the UHECR flux, but their escape through the 
CSO cocoon is made difficult by photo-disintegration. Future 
investigations of the CSO luminosity function, correlation with 
faint CSO remnants, and y-ray produced simultaneously with 
UHECRs may test the CSOs as UHECR sources. 
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